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The Chern number, a topological invariant characterising the ground state of fermionic systems, 
is a definitive theoretical signature that determines whether a given superconducting system can 
support Majorana zero modes. Here we show that it can be faithfully reproduced through a set of 
time-of-flight measurements, making it a diagnostic tool also in actual cold atom experiments. We 
demonstrate the measurement scheme by employing a chiral topological model of spinless fermions 
that only requires the experimentally accessible s-wave pairing and staggered tunnelling that mimics 
spin-orbit coupling. By adiabatically connecting this model to Kitaev's honeycomb lattice model, 
we show that it gives rise to v — dbl phases, where vortices carry Majorana fermions, and v — ±2 
phases that emerge as a collective state of localised Majorana fermions. Hence, the preparation 
of these phases and the detection of their Chern numbers provides an unambiguous signature for 
the presence of Majorana modes. Finally, we demonstrate that our detection procedure is resilient 
against most inaccuracies in experimental control parameters as well as finite temperature. 



PACS numbers: 67.85.-d,03.65.Vf 

The discovery of a new particle carries an excitement 
that reaches well beyond the scientific community. Ma- 
jorana "fermions" pQ, even though they are expected to 
emerge as zero mode quasiparticles in topological sys- 
tems, are highly anticipated. This is due to their pre- 
dicted exotic statistics that is unlike of any other be- 
haviour we have seen so far. Importantly, they offer the 
possibility of employing their statistics for sophisticated 
technological applications such as topological quantum 
computation [2]. 

For a fermionic system to support localised Majorana 
zero modes, two conditions need to be satisfied: the 
spectrum should possess particle-hole symmetry and the 
ground state should exhibit suitable topologically non- 
trivial behaviour. Particle-hole symmetry implies that 
for a stationary state ^ E with energy E, there exists an- 
other state ty-E with energy — E. The suitable topolog- 
ical character of the ground state necessitates the pres- 
ence of chiral edge states [3], which in turn imply that 
zero energy modes can be localised at the core of vortices 
[4]. Due to particle-hole symmetry these E = modes 
satisfy ^\ — i-e. they describe excitations that are 
their own anti-particles, known as Majorana zero modes. 

Particle-hole symmetry is a natural property of super- 
conducting fermionic systems. To detect whether a given 
system can support Majorana zero modes is therefore 
equivalent to detecting the topologically non-trivial na- 
ture of the ground state. This can be given in terms of a 
topological quantum number v, the Chern number, that 
takes integer values [5 . All systems with non-zero v have 
chiral edge states, but only systems characterised by odd 
v support localised Majorana modes [6 . However, the 
detection of the Chern number is challenging due to its 



non-local nature. Examples exist where this topological 
invariant is directly related to physical observables, such 
as the off-diagonal conductivity in the quantum Hall ef- 
fect [7] , but in other microscopically distinct systems such 
connections are in general unknown. 

In this work we show how to robustly measure the 
Chern number of a superfluid of fermionic atoms from 
time-of-flight images [8, 9]. We demonstrate our detec- 
tion method in the context of an experimentally tractable 
superconducting model that is adiabatically connected 
to Kitaev's honeycomb lattice [6]. The model supports 
phases with Chern numbers v = 0, ±1 and ±2. The 
v = ±2 phase emerges uniquely as the collective state of 
a lattice of Majorana fermions [20] , Thus, the experimen- 
tal preparation and detection of these phases constitutes 
two types of direct evidence for the existence of Majo- 
rana modes: (i) detecting v = ±1 phases identifies that 
localised Majorana modes can be supported in the system, 
while (ii) detecting the v = ±2 phases is an unambiguous 
signature that Majorana modes do exist in the system. 
We show that both the simulation of the superconduct- 
ing model and the required time-of-flight measurements 
can be robustly implemented in state-of-the-art ultracold 
atom experiments [T0HT2] . 

Chiral topological order with s-wave pairing:- The re- 
alisation of stable p-wave pairing, the usual route for Ma- 
jorana modes in superconducting systems [5], has been 
found to be challenging with cold atoms [13]. As an al- 
ternative, various routes of combining an effective spin- 
orbit coupling together with s-wave pairing have been 
proposed p~4j [15] . All these proposals rely on synthesis- 
ing the spin-orbit coupling using multiple atomic internal 
states. Here we show that qualitatively similar spin-orbit 
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FIG. 1: A chiral topological superconductor with s-wave 
pairing on a square lattice. When tunnelling of constant 
magnitude \t\ is along the arrows, the numbers denote the 
complex tunnelling phases. This pattern can be achieved by 
employing Raman assisted tunnelling with a different set of 
beams for each direction. The pairing is between neighbour- 
ing atoms in two distinct internal states, denoted here as cir- 
cles and squares, which facilitates the experimental imple- 
mentation of a real uniform pairing of magnitude A. Black 
(white) columns denote lattice sites with chemical potential 
lib — fi -\- 5 (fiw = fi — S), where the detuning 5 is controlled 
by a columnwise energy offset Voff- This overall coupling 
pattern is translationally invariant with respect to a two site 
cell containing a black and a white site. 
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FIG. 2: Left: The phase diagram of model |T]). Colour 
encodes the magnitude of the fermionic spectral gap, as a 
function of the overall chemical potential /i and the detuning 
5. The Chern number v for each phase is shown. The phase 
diagram is fully symmetric with respect to /i —> —/i, while for 
5 — »> —5 all the Chern numbers become time reversed (y — »• 
—v). As shown in Appendix D, the regions fi < 5/2 (/i > 5/2) 
can be identified with the honeycomb model in the presence 
(absence) of a vortex lattice. Right: The winding number 
v (encoded in colour), as obtained from the observables ^5j, 
shows excellent agreement with the Chern numbers v. The 
plots are for t — 1 and A = 2. 
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physics can be achieved with a single atomic state by 
suitably braking translational symmetry. 

Our model is defined for fermions on a square lattice 
and combines staggered complex hopping with a uniform 
superconducting s-wave interaction. The Hamiltonian is 
given by 
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where aj creates a fermion at site j = (j x ,jy), the tun- 
nelling amplitude t and pairing potential A are both real 
and the chemical potential /ij = \i + (— iy x 5 is stag- 
gered by the detuning 5. Translational symmetry is bro- 
ken along the x-direction with the "magnetic" unit cell 
(there is no net flux) consisting of a black and a white 
site, as shown in Fig. [T] Inspired by the Kogut-Susskind 
fermions [I6j Hi] we interpret this lattice degree of free- 
dom as a "pseudospin" r G {6, w} of the fermions a\ y 
As detailed in Appendix A, we can thus view the Hamil- 
tonian as an effective pseudospin-orbit coupled sys- 
tem: Tunnelling along the ^(y)- direct ion changes (con- 
serves) the pseudospin state, which effectively realises an 
anisotropic Rashba type spin-orbit coupling, while the 
chemical potential detuning plays the role of a Zeeman 
term. Thus, by adding s-wave pairing, one expects to 
find Majorana mode supporting topological phases [26] . 
To verify this, we solve ([I]) by Fourier trans- 



H= / ^ff(p)V p d 2 p, 
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where the Brillouin zone (BZ) spans p x G [0, n] and 
p y G [0, 2tt] and the Bloch Hamiltonian H (p) is a 4 x 4 
matrix. The analytic solution to the model is presented 
in Appendix A. Fig. [2] depicts the phase diagram, which 
shows that by varying the chemical potentials we can 
move between a variety of chiral topological phases with 
Chern number v = 0, ±1 and ±2. Each phase is sup- 
ported by a large range of parameters, thus facilitating 
its experimental realisation. 

As we explicitly show in Appendix B, our model is adi- 
abatically connected to Kitaev's honeycomb model [6], 
where chemical potential staggering can be related to 
the presence of a background vortex lattice. This con- 
nection implies that in the non-Abelian v = ±1 phases 
(in the weakly staggered /i > 5/2 regimes) isolated vor- 
tices will bind localised Majorana modes with short range 
interactions [27]. However, when these vortices form a 
regular array (the strongly staggered \i < 5/2 regime 
[19]), the interactions cause the Majorana modes to hy- 
bridise, which has been shown to lead to a unique collec- 
tive Abelian state with v = ±2 [20]. Therefore, by sim- 
ply tuning the chemical potential detuning we can move 
between regions with and without underlying vortices. 
Observing the corresponding changes in Chern number 
provides thus direct evidence for the presence of Majo- 
rana modes through their collective behaviour. 

Optical lattice implementation:- To physically realise 
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([T]) we suggest to use two interpenetrated 2D optical 
square lattices. Each lattice hosts a different internal 
state of the same species of neutral fermionic atoms 
denoted in Fig. [I] by circles and squares. The lat- 
tice is deep enough to suppress ordinary hopping along 
all three spatial directions. Confinement along the z- 
direction is taken to be even deeper, so that each xy- 
plane can be regarded as an independent copy of the 
same experiment. A Feshbach resonance activates the s- 
wave pairing p~4j [21] between atoms in different internal 
states, using the fact that the spacing between sublat- 
tices is smaller than the separation between atoms in 
the same state. Hopping between sublattices is achieved 
by Raman-assisted tunnelling [22-24 , a technique that 
allows both for complex tunnelling amplitudes and also 
for a delicate control of the chemical potentials [see Ap- 
pendix C]. The hopping parameter structure can be im- 
plemented by means of an energy offset, Voff, hi alter- 
nate columns. Note that the complexity introduced by 
having two species copes with the fact that it is diffi- 
cult to achieve a Feshbach induced s-wave pairing with 
atoms in the same state. Should the experiment allow 
for a same-species Feshbach coupling, the whole setup 
could be simplified using just a single lattice and the flux 
inducing mechanism from Ref. [25] . 

Observables for the Chern number:- We now show how 
to construct a quantity that (i) is an integer, (ii) is defined 
in terms of measurable quantities and (iii) reproduces the 
Chern number in the zero temperature limit. Formally, 
the Chern number, is defined as the winding number 
of the projector onto the ground state [7]. When the 
Bloch Hamiltonian is a 2 x 2 matrix, i.e. the system is 
fully translationally invariant, it can always be written 
as H(p) oc S(p) • cr for some vector field S(p). Here a 
denotes a vector of Pauli matrices. The Chern number, 

is then equivalent to the winding number 

4t<" Jbz |S(p)| 3 V d Px dp y J 

(3) 

which counts how many times S(p) winds around the 
Bloch sphere in the particle-hole space as one spans the 
whole Brillouin zone [7 . We can obtain this quantity 
from the components of the vector field S(p). These 
components are observables that can be obtained as the 
ground state expectation values 

S(p) = <*| S p I*) , S p = V p <tV p , (4) 

with the physical observables X) p being given in the basis 
i/?p = (a p ,a_ p ) of the Bloch Hamiltonian H(p). 

To generalize the observable Q to staggered sys- 
tems with m site unit cells, we define an indepen- 
dent vector field S^)(p) = (^| p |^) for each 
of the i = l,...,m sublattices. When the full 
2m x 2m Bloch Hamiltonian is given in the basis 



= ( a i, P > ■ ■ ■ > a L,P' ai -P' ' • • ' a m-p) , the correspond- 
ing sublattice observables are 5](i), p = p "*!*^)^ w ^h 

P = ( a (i) P ' a (i)-v)' Explicitly, their components are 
given by 

S fi),p = a (i),p a (i),-p + a (i) -P a (i),p> 

S (i),p = -H),p4),-p +ia (*)'-p a ( i )'P' ( 5 ) 

S (i),p = a (i),p a (i),P - a (i),-p a (i),-p- 

Substituting these into Q, we can construct m winding 
numbers v<^ — £[S^)], i = l,...,m. We now define the 
total winding number by 

m 

By direct comparison we observe that, for the case of the 
two site staggered model , the winding number v is in 
excellent agreement with the Chern number v for all cou- 
pling regimes, as shown in Fig. [2] In particular, phases 
with different \v\ are exactly distinguished, so the observ- 
ables ^b,p and T> w p can be used to faithfully identify the 
topological phases of our model [see Appendix D]. 

The experimental measurement of 5](i), p proceeds in 
several steps. Time-of- flight images give direct access 
to (aj^ p a (i),p)i f° r both i — b,w separately, which is 
enough to obtain S z . To obtain the orthogonal com- 
ponents S x and S y we note that pairing in the Hamil- 
tonian ([TJ along the ^/-direction is proportional to the 
operator T, y . An evolution under only this term imple- 
ments a rotation around the S y axis, which maps the 
value S x onto the measurable S z . Thus, to measure S x , 
we suggest first applying this rotation for a brief period of 
time, for instance by ramping up the optical lattice along 
the x-direction and switching off the hopping, and only 
then performing the time-of- flight measurement. While 
the setup contains no similar rotation to map S y to S z , 
it can be obtained from noise correlation measurements 
once both S x and S z are known [see Appendix D]. As 
the time of flight images are digitalised one could either 
compute the discretised version of ([3| [9 or perform a 
continues fit to the values of S x , S y and S z . 

Robustness to experimental errors:- Our measurement 
procedure is robust under the two main cold atom ex- 
periment imperfections - a harmonic trapping poten- 
tial and a finite temperature. We simulate the first 
in a finite L x L system with open boundary condi- 
tions by introducing in the chemical potential /ij = 
fji + (-1)1*5 + Md 2 cu 2 [(j x - L/2f + (j y - L/2) 2 ], where 
M is the mass of the atomic species and d is the lattice 
spacing. Such a potential removes the idealised condi- 
tions of translational symmetry that is required for the 
formal definition of the Chern number. Assuming that a 
local density approximation holds, this spatially depen- 
dent chemical potential induces the coexistence of dif- 
ferent phases: some of insulating character, some not; 
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FIG. 3: The effect of a harmonic trap for the atoms on 
the winding number v. Left: Winding number v (encoded 
in colour) as measured by the observables Qinal5xl5 
lattice with a trapping frequency fvuj/t = 0.5. Right: The 
winding numbers as the trapping frequency uj is increased for 
the v = 2 phase (<5,/i)=(3,0), [black solid line], the v — 1 
phase (£,/i)=(4,2), [red dashed line] and the z/ = — 1 phase 
(J,/x)=(l,3), [blue dotted line]. Both figures assume a lattice 
spacing of 500nm and the atomic mass of 40 K. 

some with topological order, some with no order at all. 
Nevertheless, our time-of-flight images will suppress the 
contribution from all trivial phases. 

Comparison between Fig.[2]and Fig.[3]shows that while 
the phase diagram is deformed, all the topological phases, 
as characterised by the winding number z>, remain robust 
for an extended range of trapping frequencies uj. Since 
the same winding number reproduced the Chern numbers 
in the idealised unperturbed system, we regard that it 
faithfully capturing the nature of the topological order 
also in perturbed systems. 

To model the effect of finite temperature T we restrict 
to fermionic excitations in the lower band with no ther- 
mal vortex excitations. Then the ground state |^(p)) is 
replaced by a thermal state, which has occupation n_ in 
the ground state and occupation n + in the first excited 
state. Here n± = l/(exp(± ^T^ ) + 1) * s the Fermi dis- 
tribution for the lowest lying energy band 22i(p). Eval- 
uating the expectation values Q on this state we find 
S(£) — >• S(i)T = ( n - ~ n +)S(i)- This implies that finite 
temperature only changes the normalisation of the vector 
fields. Thus, the winding number ([3| still faithfully repro- 
duces the Chern number. The decrease in the visibility 
|S^) T | 2 means that the time-of-flight measurements be- 
come more susceptible to noise, but this can be compen- 
sated by sufficient experimental precision. This contrasts 
with the behaviour of topological entanglement entropy, 
which provides an alternative probe for topological or- 
der in cold atom systems [28]. In the thermodynamic 
limit topological entropy vanishes at any finite tempera- 
ture rendering its applicability challenging [29]. Further- 
more, unlike our method, topological entropy can not 
distinguish topological phases with same total quantum 
dimensions [2]. 

Conclusions:- We have presented a general method 
to detect the Chern number of superconducting mod- 



els from time-of-flight images. We showed that our 
method can detect phases with Chern numbers v — ±1 
or v — ±2, thereby identifying those that can support lo- 
calised Majorana zero modes or arise through their collec- 
tive behaviour, respectively. Both types of phases arise 
in the proposed model, which establishes a new route 
to topological phases with cold atoms. Staggered tun- 
nellings of spinless fermions realise an effective spin-orbit 
coupling, which, in analogy to topological semiconductor 
heterostructures [26], can give rise to Majorana modes 
when s-wave pairing is introduced. We showed this ex- 
plicitly by adiabatically connecting the staggered model 
to Kitaev's honeycomb model [6], where both the proper- 
ties of the non-Abelian vortices in the v = ±1 phases [27] 
and their collective behaviour [20 has been extensively 
studied. While we demonstrated our experimentally ro- 
bust measurement scheme in the context of a particular 
model, our method is readily applicable also to other cold 
atom realisations of topological phases [I4j [15] . An im- 
portant open question is the interplay between different 
phases in an experiment with inhomogeneities, where the 
boundaries are not so clearly defined. 

Acknowledgements:- JKP would like to thank Wolf- 
gang Ketterle for inspiring conversations. This work was 
supported by EPSRC and by Spanish MICINN Project 
FIS2009-10061, FPU grant No.AP 2009-1761, CAM re- 
search consortium QUITEMAD S2009-ESP-1594. 



[I] E. Majorana, Nuovo Cimento, 5, 171 (1937). 

[2] J. K. Pachos, Introduction to Topological Quantum Com- 
putation, Cambridge University Press (2012). 

[3] Y. Hatsugai, Phys. Rev. Lett. 71, 3697 (1993). 

[4] C. Chamon et a/., Phys. Rev. B 81, 224515 (2010). 

[5] N. Read and D. Green, Phys. Rev. B 61, 10267 (2000). 

[6] A. Kitaev, Ann. Phys. 321, 2 (2006). 

[7] D. J. Thouless, et a/., Phys. Rev. Lett. 49, 405 (1982). 

[8] E. Zhao, N. Bray-Ali, C. J. Williams, I. B. Spielman and 
I. I. Satija, Phys. Rev. A 84, 063629 (2011). 

[9] E. Alba et al, Phys. Rev. Lett. 107, 235301 (2011). 

[10] E. Altman, E. Dernier and M. D. Lukin, Phys. Rev. A 
70, 013603 (2004). 

[II] S. Fo'lling et al, Nature 434, 481 (2005) 

[12] M. Greiner, et a/., Phys. Rev. Lett. 94, 110401 (2005). 
[13] J. P Gaebler, et al, Phys. Rev. Lett. 98, 200403 (2007). 
[14] L. Jiang et a/., Phys. Rev. Lett. 106, 22 (2011). 
[15] C. Zhang, et a/., Phys. Rev. Lett. 101, 160401 (2008); 

C. Zhang, Phys. Rev. A, 82, 021607 (2010); S. Zhu et al, 

Phys. Rev. Lett. 106, 100404 (2011). 
[16] P. Maraner and J. K. Pachos, Phys. Lett. A 373, 2542 

(2009) . 

[17] J. B. Kogut, and L. Susskind, Phys. Rev. D 11, 395 

(1975); L. Susskind, Phys. Rev. D 16, 3031 (1977). 
[18] D. A. Ivanov, Phys. Rev. Lett. 86, 268 (2001). 
[19] V. Lahtinen and J.K. Pachos, Phys. Rev. B 81, 245132 

(2010) ; J. K. Pachos, Ann. Phys. 322, 1254 (2007). 
[20] V. Lahtinen, et a/., Phys. Rev. B 86, 075115 (2012). 
[21] G. M. Bruun and C. J. Pethick, Phys. Rev. Lett. 92, 14 



5 



(2004). 

[22] D. Jaksch and P. Zoller, New J. Phys. 5, 56 (2003). 
[23] F. Gerbier and J. Dalibard, New J. Phys. 12, 033007 

(2010) . 

[24] L. Mazza et al, New J. Phys. 14, 015007 (2012). 

[25] M. Aidelsburger et a/., Phys. Rev. Lett. 107, 255301 

(2011) . 

[26] J. Alicea, Phys. Rev. B 81, 125318 (2010). 
[27] V. Lahtinen, New. J. Phys. 13, 075009 (2011). 
[28] DA. Abanin and E. Dernier, Phys. Rev. Lett. 109, 
020504 (2012). 

[29] C. Castelnovo and C. Chamon, Phys. Rev. B 76, 184442 
(2007); S. Iblisdir, D. Perez-Garcia, M. Aguado and J.K. 
Pachos, Phys. Rev. B 79, 134303 (2009). 



Appendix A: Chiral topological order with s-wave 
pairing 

In this Appendix we first give the analytic solution to 
our model ([I]) and then discuss its interpretation as an 
effective spin-orbit coupled system with s-wave pairing. 
We also show how sign staggered couplings in a p-wave 
model allow for the complex pairing to be traded for real 
5-wave pairing with complex hopping. 



Analytic solution 

The Hamiltonian ([I]) can be Fourier transformed with 
respect to the magnetic unit cell. Writing it in the 
particle- hole basis ip P = (a b , ttJy jP , a& 5 - p , a Wj —p), we ob- 
tain the quadratic Bogoliubov-de Gennes Hamiltonian 

h = Ibz ^i H (p)^p d<2 p^ where 



with 



H(p) = 



f± 



( /+ W*+ ih 9- \ 

-ig+ f- ~9- ih 

-ih -g*_ -/+ ig*_ 

\ 9- -ih -ig+ -/_ J 



(/i ± S) + 2tcos(p y ), 



(7) 
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h = 2Asin(p 2/ ). 
By employing a Bogoliubov transformation 

iE(b,w) 

for some complex coefficients u n i and v n ^ the Hamilto- 
nian can be brought to the diagonal form 



E n (p) b Lp b n,P + E n (p)b n - p b, 



n=l,2 J BZ 

The four energy bands are given by 



n,— P u n, — p 



Et(p) = ±\M(P) + (-1)V^ 2 (P)-4B(P), (8) 






FIG. 4: The spectral flow on a cylinder for the v = (/x = 
0,£ = 5), v — 1 (n = 2,6 — 5) and v — 2 (fi = 0,5 — 2) phases 
showing edge states crossing zero energy per edge. The 
plots are for t = A = 1. 
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The particle-hole symmetry present in the problem is 
represented by C = a x ® 11 that swaps the creation and 
annihilation operators of opposite momenta. It acts on 
Q as CH(p)C~ 1 = -iJ*(-p), which implies that zero 
energy eigenstates at the momenta p = (0,0), (0,7r) will 
be self-conjugate. Fig. [4] shows that the edge states in- 
deed cross zero energy at these momenta implying that 
they are (dispersing) Majorana modes. In the presence of 
a vortex (a puncture in the plane with 7r-flux through it), 
the edge states will become localised Majorana modes [4]. 
Odd number of edge states implies that an isolated Ma- 
jorana mode will always remain localised at zero energy 
at the vortex core, while an even number of them leads 
to complete hybridisation with all the Majorana modes 
pairing up to localised Dirac fermions. 



Staggered tunnelling as an effective spin-orbit 
coupling 

Even if there is no net magnetic field through the sys- 
tem, staggering in the couplings leads to an enlarged unit 
cell. For the coupling pattern shown in Fig. [I] the unit 
cell consists of two sites, which we denote as black (b) and 
white (w). Interpreting this lattice degree of freedom as 
a pseudospin r G w) of the fermions aj r , using the 
spinor" jbj = (a&jjG^j) we can rewrite the different 



terms of (ll| in the following way: 



(9) 



where the Pauli matrices r a act on the psedospin degree 
of freedom. This suggests the following interpretation in 
terms of the fermions iJ;J: 
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• \± still acts as the chemical potential, while the de- 
tuning acts now effectively as a Zeeman term of 
magnitude V z = S. 

• Tunnelling in x-direction realizes an anisotropic 
Rashba type spin-orbit coupling p x cr y of magnitude 
a = t and a transverse magnetic field of magnitude 

Vy = t. 

• Pairing will be of uniform amplitude A, but it will 

be an anisotropic mixture of singlet (x-direction) 
and triplet pairing (^/-direction). 

The elements above - the spin-orbit coupling, magnetic 
fields of different direction and the s-wave pairing - are 
the components of Majorana mode hosting semiconduc- 
tor heterostructures [26] . 

These components appear also in a staggered p-wave 
superconductor, where, similar to our model, the stag- 
gering allows complex pairing to be traded for complex 
hopping. Denoting the fermionic operators aj by frj (wj) 
when they live on the black (white) sublattice, shown in 
Fig. [I] we can write the simplest possible p-wave lattice 
model as [1 



TTU 

Px+iPi 



= J2 [(ji+b]bi + (i-wjwj) 



( 



+t ( ub\b 



'j+y +W\Wj+y 



ib\w] 



iw]b 



'j+x 



(10) 



-ub\b\ +9 



w l w l+y 



b\ 



■4 h U 



where fi± = [i ± 5, t and A are again the tunnelling 
amplitude and the pairing potential, respectively. To 
stagger the couplings, we have introduced the parame- 
ter u. While u = 1 recovers the original Hamiltonian of 
PQ (up to a unitary frj — >• zfrj and subsequently identifying 
bj = Wj) : setting u = — 1 breaks translational symmetry 
by flipping the signs of the tunnelling and pairing on the 
black sites. This results in the pairing becoming real and 
uniform s- wave pairing, while the complex p-wave sym- 
metry is transferred to the tunnelling amplitudes. In ad- 
dition, the sublattice dependent sign staggered tunnelling 
in the ^/-direction means that we have introduced 7r-flux 
on each square plaquette. Fourier transformation gives a 
Hamiltonian of the form ^ with the only difference that 
one needs to replace /+ — >• = (/i + 8) + 2utcos(p y ). 
Diagonalisation gives a phase diagram that under swap- 
ping fi <H> S is identical to that of Fig. [2j Staggering 
enables thus to trade complex pairing into complex hop- 
ping (with or without net magnetic field), while preserv- 
ing topologically non-triviality the phase diagram. This 
can be advantageous from the point of view of experi- 
ments where complex tunnelling is in general easier to 
realise than a complex pairing. 

Writing (10) in terms of the spinor ijjj = (&J,wt), we 
obtain a similar interpretation of the Hamiltonian de- 
scribing a pseudospin-orbit coupled system. The only dif- 
ference to d9j) is that the sublattice dependent tunnelling 



in ^/-direction gives rise to a further spin-orbit term of 
Dresselhaus type (ipj T z ifjj+ y ). Borrowing from the be- 
haviour of the topological heterostructures [26j, this ex- 
tra term enables us to gain insight into the swapped be- 
haviour of the models when tuning \i and S. There it is 
known that in the presence of only Rashba coupling, i.e. 
in the original model 0, a finite perpendicular Zeeman 
field V z = S is required to drive the system into a topo- 
logical phase. On the other hand, when a Dresselhaus 
coupling is present, an in-plane magnetic field V y = t is 



sufficient. Indeed, the model (10) above has topological 
phases even for 5 = 0. 



Appendix B: Adiabatic connection to Kitaev's 
honeycomb lattice model 

In this Appendix we demonstrate that our staggered 
model is adiabatically connected to Kitaev's honey- 
comb model [6]. We show this explicitly for the v = ±2 
phases that arise in the full- vortex sector of the model [19] 
due to the collective behaviour of the Majorana modes 
bound at the vortex cores [20] . 

The honeycomb model is a local spin lattice model that 
can be adiabatically connected to the p x + ip y super- 
conductor [2 . It can be solved exactly by mapping it 
to a tight-binding model of free Majorana fermions on 
the honeycomb lattice with nearest (magnitudes J x -,J y 
and J z depending on link orientations) and next nearest 
neighbour (magnitude K) tunnelling [6] . After this map- 
ping, the Hamiltonian in the full- vortex sector (7r-flux on 
each hexagonal plaquette) can be written as [2 [19] 

j 

+ iK ^2(-l) Jx [ajttj-x + OjOj+y + frjfrj+x + &j&j-y] , 



H.c. 



(11) 



j cluv^ ^ — bj are Majorana operators on 



where a\ = a\ and b\ = 
j j j 

the two triangular sublattices of the honeycomb lattice. 
The simplify the model we have included only four out 
of the six possible next nearest neighbour hoppings (the 
omitted terms being ajOj+x+y and 6j6j +x+y ). 

The full- vortex sector is encoded in the (— l) Ja = factors 
that stagger the signs of the Majorana hopping ampli- 
tudes J z and K. Pairing the Majorana operators into 
complex fermions Cj by 



i 



(12) 



(the phase 0j to be defined below) the links connecting 
the Oj and bj sites of the honeycomb lattice become the 
sites of a square lattice. In this basis the Hamiltonian 
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FIG. 5: When the Major ana operators Oj and frj are paired 
into complex fermions Cj, the links connecting them become 
sites of the square lattice. The resulting Hamiltonian with 
couplings ( 14 ) has a four site unit cell shown here as the 



dashed boxes. Redefining the operators on the circle (square) 



sites by Cj 



j( c J 



^Cj) preserves the real pairing 



for arbitrary x while unitarily transforming the hopping am- 
plitudes. For x = 0/2 — 7r/4 they are brought to the form 
(15). The chemical potential /ij alternates in the x-direction 
between 2J Z (black sites) and — 2J Z (white sites). 



takes the BCS form 
j 

A ^CjC j+ x + A^cjCj+y] + H.c, (13) 



where we defined 



2J z (-l) j * 



t x = re i(-iy*(2-(-iyy) 
t y _ 2Je~ i (- i y x+jy( t ) 



(14) 



A x = 2J, 



r, 



with J = J x - 
arctan(J/(2iT)). 
12k is given by #j 



m 



J y , r = v/(2J) 2 + (4if) 2 and = 
In terms of these variables the phase 



le variables 



_(_l)j, 
r and (7> 



are all functions of the 



honeycomb couplings J, J z and if. From now on we will 



treat them as independent variables and show that ( 13 ) 
can be adiabatically connected to ([I]). First of all, we set 
the scale by directly identify the chemical potential de- 
tuning to be S = 2 J z with the overall chemical potential 
set to /i = 0. Next, we would like to tune r — > 2 J = £ 
in order to make the tunnelling and pairing amplitudes 
equal. Fig. [6] shows that the gap increasing during this 
process. Finally, we would like to tune the phase of the 
tunnelling amplitudes (14). As illustrated in Fig. |5j these 
couplings are translationally invariant with respect to a 




(1-xH + ti/4 x 



0.2 0.4 0.6 0.8 10 0.2 0.4 0.6 0.8 1 

X X 

FIG. 6: Adiabatic connection between the v — — 2 phases 
in the full- vortex sector of the honeycomb model and in |l]). 
Left: As we first tune r —t 2 J for the tunnelling couplings 
t* ,y and the pairing potentials A a to have equal amplitudes, 
the gap A increases. Right: Tuning subsequently cj) — >> ty/4 to 
match the phase pattern in the hoppings, we find the gap first 
increasing and then settling to a constant value. The data is 
for J — J z — 1 and K — 0.1. Both transitions are performed 
with a linear ramp parameterised by x £ [0, 1]. 



four site unit cell and there are altogether six indepen- 
dent tunnelling couplings for a general cj). These can be 
unitarily transformed to the form 



4-X 
l l 
4-X 

l 2 

4-X 

r 3 

4-X 

l 4 



l 2 



te -in/2 
£g4i</>— i7r/2 
£ e -4i</>+i7r/2 
^ e -2i0+i7r/2 
te 2icf)-i7r/2 



it, 
-it, 

— )► it, 
-it, 

-+ t, 



(15) 



and subsequently tuned to the special value (j) = 7r/4 
(the — » above), which gives a two site unit cell with 
the same tunnelling phases as in the staggered model 
0. Fig. [6] shows the gap remains non-zero during 
this process, which implies that the models are adia- 
batically connected. Indeed, for J = J z = 1 and 
K > the honeycomb model is in a nucleated v = — 2 
phase [19]. These honeycomb couplings correspond to 
t = A x = A y = —(5 = 2 and /i = for which, as shown 
in Fig. [2j we also find a v = — 2 (the chemical potential 
detuning 5 is here negative). 

A remarks is in order. Since we put the sign staggering 



in (11) on the couplings J z and subsequently paired the 
Major anas along these links, the resulting square lattice 
model has a staggered chemical potential, but no flux 
per plaquette (one can verify that the phases in the tun- 
nellings (14) cancel). This is basis dependent though. 
Had we paired the Majoranas along, say, J y couplings, 
we would have arrived at a model with a uniform chem- 
ical potential and 7r-flux per plaquette. This shows that 
there is nothing special about having a real magnetic 
field in the system - it is the staggering that encodes the 
physics. 

The full phase diagram of the full- vortex sector of the 
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□ o 

UxJy) UxJy + l) 

FIG. 7: Energy level structure, showing how two sets of Ra- 
man beams allow for independent transitions due to the pres- 
ence of the offset Voff- This offset can be tuned to allow 
for different values of the difference in chemical potentials 
fiw — fib- The indices (j x ,jy) stand for horizontal and vertical 
position on the lattice, respectively. 

honeycomb model has been studied in [19 . We can em- 
ploy it to understand some of the features of the phase 
diagram of Fig. [2j First of all, in the dimerised limits 
\J Z \ ^> | J | one always finds a v — phase, in agree- 
ment with us finding v = phases in both 5 ^> fi and 
/i ^> 5 limits. For intermediate values one is either in a 
1/ = ±1 or a ^ = ±2 depending on the degree of stag- 
gering 5 = 2J Z . While the honeycomb model has always 
H = chemical potential when derived form the spin 
model, we can introduce a finite \i by hand. Doing so 
one finds that for \i > 5/2 the non-Abelian v = ±1 phase 
is recovered even if some sign staggering remains [19]. 
This is in agreement with Fig. [2j which shows along the 
/i w 5/2 line a similar transition between the Abelian 
v = 2 and the non-Abelian v — — 1 phases. We thus 
interpret the region /i < 5/2 as the one where the col- 
lective vortex lattice effects, and those of the Majorana 
modes bound to them, are relevant. The \i > 5/2 re- 
gion, however, should be understood as having similar 
physics as the vortex-free sector where isolated vortices 
bind localised Majorana modes. 




FIG. 8: Illustration of the vector field S for a v — — 1 
phase. Here we plot the values of S(&)(left) and S^) (right) 
for (5, /x)=(l,3). It can be seen that S(&) winds once around 
the origin, thus giving a partial \i>b\ = +1 contribution, while 
S(^) does not enclose the origin so it gives zero contribution. 

lows for independent Raman tunnelling in each direction 
and acts as a knob for changing the value of the chemical 
potential difference fi w — fib, as it can be seen in Fig. [7] 
Our implementation employs a Raman- assisted tun- 
nelling scheme on an optical lattice with a pattern of 
phases, as shown in Fig. [I] 221(24], In this figure we 
use a convention that the tunnelling direction is set to 
go from one species (circle) to the other (square). Re- 
versing the direction complex conjugates the tunnelling 
amplitude. The tunnelling element between sites j and j' 
in a Raman transition assisted by two lasers of wavevec- 
tors ki and k2 of amplitude V$ can be parameterised 
as t = Voe lcl ' r+ to(d, q), where to is a real number which 
only depends on the nearest-neighbour distance d and 
the difference between the Raman beam wavevectors 
q = ki — k2. Also r + = (j + j')/2 is the midpoint 
between the two neighbouring sites. The phase of the 
hopping parameter is thus determined by wavevector q. 
We can see in Fig. [I] that two different phase wavelengths 
for the horizontal and vertical transition amplitudes are 
needed. So implementation of this hopping pattern re- 
quires a Raman pair for each direction. Fig. [7] shows the 
energy level structure which accomplishes this pattern. 



Appendix C: Optical lattice implementation 

We propose to generate the lattice of model ([!]) by fo- 
cusing the diffracted image from a holographic mask onto 
the focal plane of an extremely large aperture lens [4]. 
The sublattices host different hyperfine states of the same 
atomic species which are sensitive to different polarisa- 
tions of the trapping laser beams. Both sublattices can 
be displaced one on top of another by means of an electro- 
optic modulator [5]. 

We also superimpose a superlattice potential V = 
Voff sin 2 (7nr /2d) whose amplitude can be controlled by 
changing the intensity of the laser beam, which cre- 
ates this potential. This superlattice structure effectively 
adds a tunable energy offset Voff to every other column 
in the lattice. This offset serves a two- fold purpose: it al- 



Appendix D: Measurement protocol 

To obtain the winding number v one needs to first con- 
struct the two independent integer- valued winding num- 
bers z>5 and v w . The structure of the observable vec- 
tor fields S(5 ?w ) itself also provides valuable information 
regarding the topology of the system. Fig. [8] shows a 
3D-representation of these vector fields in a state with 
v = — 1. We can see that a non-trivial contribution ap- 
pears whenever the vector field structure winds around 
the origin (it is assumed that the vector field is nowhere 
zero, which is reasonable if the temperature is kept low 
enough). These diagrams may hint a non-trivial topol- 
ogy of the system even if the winding number cannot be 
accurately computed due to experimental imperfections. 

Our characterisation of the quantum phases is based 
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on the winding of the vector fields S^), which are de- 
fined independently for the "black" (i = b) and "white" 
(i = w) sublattices distinguished by their different chem- 
ical potential. The energy offsets between them can be 
employed to release atoms from one of the two sublat- 
tices at a time, thus obtaining n^p and n WjP , and the 
observable Sfo(p) = ^(i),p + ^(i),-p — 1 f° r both i = b,w. 

For the components S x and S y we suggest employ- 
ing two complementary techniques. The reason for this 
is that measuring the transverse components requires 
momentum-dependent operations and having two differ- 
ent methods allows for a more robust calibration and 
redundancy in the measurements. 

The first method takes the same time-of-flight images 
which are used for S^(p) but processes them to obtain 
the noise correlations, (n^^n^ _ p ). Using Wick's theo- 
rem we can show that such an observable can be written 
in the form 

(^(i), P ^(z) -p> = ( n (i), P X n w -p> + \( a \i) P a \i) - P )l 2 

(16) 

+ ( a (i),p a W-p)( a W,p a (t)-p^ 

We can eliminate in this equation the values 
( a (i),P a \i) _ p )? f° r they can be obtained from time- 
of-flight measurements after Bragg scattering with 
momentum 2p; but in addition they are also zero in our 
setup. This means that in essence noise correlation give 
us access to the orthogonal projection of the pseudo-spin 
components, (S^) 2 = (S x ) 2 + (V) 2 . 

The second method consists on applying the pairing 
interactions for a very brief time prior to the time-of- 
flight images. The Hamiltonian when we switch off 
the hopping t = and leave only the pairing along the y 
direction, acquires the form 

H rot = Asm(p y )i (4),p4),p " Hx -) a S k P - ( 17 ) 
This operator implements a rotation around the S y axis, 



mapping the value of the S x operator onto the S z axis 

S*(p, t) = cos(0 p )S*(p, 0) + sin(# p )^(p, 0), (18) 

with 6 = Asm(p y )t/H. This allows us to extract the 
value of S x , with some difficulties around p y = 0, ±7r, to 
be discussed now. 

The value of S y is implicit in the previous measure- 
ments, through the formula, (S y ) 2 = (S^) 2 - (S x ) 2 . It 
can also be obtained experimentally using a two-step pro- 
cess. The first part consists on an evolution with hopping 
along the y direction, that maps S y to S x , while the last 
step uses again the paring to measure S x . 

In practice, as the rotation angles of the opera- 
tions mentioned above becomes infinitely slow for p y = 
0,±7r/2,±7r, we instead suggest to use a clever fitting 
procedure, extrapolating smoothly the values of S from 
the measurements of S*^ 00 ^ 1 - and S y . From our nu- 
merical experiments, this can be done effectively, as the 
norm of S does not become zero anywhere, and the an- 
gles in the xy plane behave smoothly across the Brillouin 
zone. Finally, note that in practice the field component 
S^(p) could be sufficient to distinguish between Vi = 
and Vi = ±1 by counting the number of times Sfo takes 
extremal values over the Brillouin zone. An odd num- 
ber of maxima, say three, over one minimum is obtained 
with the right configuration of Fig. [8j while even num- 
ber of maxima and minima should correspond to winding 
number zero. 
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